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■ We analyse products of random RxR matrices by means of a variant of the replica trick which was 

' recently introduced for one-dimensional disordered Ising models. The replicated transfer matrix can 

be block-diagonalized with help of irreducible representations of the permutation group. We show 
that the free energy (or the Lyapunov exponent) of the product corresponds to the replica symmetric 
jj^ | representation, whereas non-trivial representations correspond to certain correlation functions. 
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The asymptotic properties of products of random matrices play an important role in many physical problems . 
In models like disordered one-dimensional magnetic systems they describe the thermodynamic quantities such as 
free energy or correlations, for localization of electronic waves in random potentials they are related to the transport 
Ch ' properties, see also Q. Such products also appear in the context of chaotic dynamical systems characterizing the 
1 I divergence of neighboring trajectories. 

Although there are many known results on products of random matrices, some of them even mathematically 
rigorous, we want to present a general replica transfer matrix method. Replicas are known to be a very powerful but 
nevertheless somewhat mysterious tool in the statistical mechanics of disordered systems and related problems Q . In 
the case of mean-field models replicas predict the concept of replica symmetry breaking which is related to a highly 
nontrivial ultrametric structure of states in the low temperature phase. 

The existence of replica symmetry breaking in low dimensional systems is not yet clear, see e.g. the argumentation 
of H|. In H a replica approach to one-dimensional disordered Ising models was presented. Although there does not 
exist any phase transition at nonzero temperature, a rich replica structure could be observed leading to a 'natural' 
critrion for replica symmetry breaking in this special system which is not related to Parisi's replica symmetry breaking 
Q\ , scheme for mean-field models. This criterion is based on the representation structure of the permutation group and 
could be deduced to a large extent with rigorous methods. In the present paper we generalize the approach of || 
to infinite products of random matrices of finite dimension. We mainly use the language of statistical mechanics, 
i.e. the random matrices are considered as transfer matrices of one-dimensional models with finite discrete degrees 
of freedom and random short-range interactions. We show that the representation theoretic approach to replica 
symmetry breaking is quite general and can be formulated without specifying a particular one-dimensional model. 
The outline of the paper is the following. In Sec. II we introduce the replicated transfer matrix. For the analysis we 
. need several tools from the representation theory of the symmetric group. These are presented in the third section. 
In Sec. IV the replica symmetric representation space is considered and the free energy is calculated. The connection 
between non-trivial representations of the symmetric group and connected correlation functions is analyzed in Sec. 
O ■ V. They will provide a natural criterion for replica symmetry breaking. In the last section we give a summary and 
\ outlook. Several appendices contain longer calculations or proofs. 
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II. THE REPLICATED TRANSFER MATRIX 

We consider N R x R matrices Ti,i — l,...,n, drawn from a single probability distribution P(T), where R is any 
positive integer. In the case of an one-dimensional model with random Hamiltonian H = Hi(si, Sj+i) (sj can take 
R different values) and inverse temperature they are given by Tj = (exp{iJ;(si, Sj+i)}). For a general distribution 
these matrices do not commute. Therefore we cannot find a common system of eigenvectors. In order to calculate 
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self-averaging quantities such as the free energy we introduce as usual the n-fold replicated and disorder averaged 
partition function, 



= < (trJjT,)" > 



i=l 



= (tr « T®" ») iv , (1) 

where <C • 3> denotes the average with respect to P{T) and (g) the Kronecker product of matrices. With this relation 
we are able to replace the product of N random R x R matrices by the A^-th power of a single R n x R n matrix 
which can be analyzed using standard transfer matrix techniques: we have to find expressions for the eigenvalues of 
T n :=<C T® n which enable an analytic continuation in n. The free energy is then given by 

/ = -i « InZ > = -~ lim d n « Z n » , (2) 

fj fj n^O 

which is dominated by the largest eigenvalue of T n for n — > 0. Several correlation length can be described by smaller 
eigenvalues of the same matrix. 

To calculate this, some notations will be introduced. The original matrices Ti act on a i?-dimensional vector space 
V. As a basis we chose any orthonormalized set of R vectors and denote these by \s), s = l,...,R. Consequently 
T n is a linear operator defined on the n-fold tensor product V® n of V with itself which has dimension R n . The 
orthonormalized basis- vectors of this space are chosen naturally as Is 1 ) ® |s 2 } <S> ■■■ <8 \s n ) =: \s 1 s 2 ...s n ) where s a G 
{1, R} for all a = 1, n. The matrix elements of T n are then given by 

n 

(s 1 s 2 ...s n \T n \s' 1 s' 2 ...s' n ) = « n( sa l T l s ' a ) » 

a=l 
n 

= « n r SaK » (3) 

a=l 

for any two basis vectors of V® n . 

The average over P(T) produces interactions between the replicas. Nevertheless the replicas are completely equiv- 
alent, a renumbering does not change the matrix T n . This leads to a symmetry of the transfer matrix under replica 
permutations, i.e. to replica symmetry of T n . The action of any permutations is given by the R n dimensional 
representation D of the symmetric group S n : 

^Ml^s 2 ...*"} = i^VPL.W")), VnESn, (4) 

whose operator product with T n commutes, 

D(tt) T n - T n D(tt) , Vtt e S n . (5) 

A direct consequence of equation (||) is the closure of any eigenspace of T n under permutations, these eigenspaces 
define a subrepresentations of D which in the most general case are irreducible. Further reducibilities would be a hint 
to a further hidden symmetry. 

Consider an element Y of the group algebra s„ of S ni i.e. Y is a linear combination of permutations i 6 5„. 
Due to (ra) and the linearity of the action of the transfer matrix on V® n it also commutes with T n . The space 
U = Y V® n = J^s 1 S " RY\s 1 s 2 ...s n ) is therefore invariant under action of T n . If we are able to construct elements 
of s„ projecting V® n to a proper subspace we can thus achieve a block diagonalization of T n by its restriction to U 
and to its orthogonal complement (1 — Y) V® n . 



III. SOME REMARKS ON THE SYMMETRIC GROUP 



In this section we review some properties of the symmetric group and its irreducible representations. These are 
well-studied and numerous excellent presentations can be found, e.g. in 0^]. Here we omit any proofs. 

The symmetric group S n contains the n! permutations of n distinguishable objects. Consider any representation 
D on a linear space V. D is called to be irreducible iff there are no proper subspaces {0}) of V closed under 
D(S n ). A representation is called completely reducible iff it can be decomposed into a direct sum of irreducible 
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subrepresentations. This decomposition is unique up to isomorphisms. Our D defined in the previous section is 
completely reducible. 

The irreducible representations of S n are classified by the so-called standard Young tableaus. Each Young tableau 
is characterized by a partition of n, i.e. a set of integers Ai > A 2 > ... > A m > 0, m < n, fulfilling J2 a ^ a = n - One 
arranges m rows of length Ai, A TO as shown in the figure and fills the boxes with the integers 1, n. The tableau 
is called standard iff the entries of the boxes are increasing within every row and within every column, see e.g. the 
figure. 
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FIG. 1. Examples for standard Young tableaus for n — 9 



At first we define the row symmetrizer SYM[ Al ....,A m ] = EH^Li SYM a with SYM a being the sum of all permutations 

within the a-th row. Then we still need the column antisymmetrizer ASYM[ Al Xm ] = Ilbii ASYMf, with ASYMb 

being the total antisymmetrizer of column 6, i.e. the sum of (— I) 71 "^ over all permutations of this column. (— I) 71 " 
signifies whether tt is odd or even. The Young operator is then defined by 

Y [Xu ..., Xm] = ASYM [Ali ..., Am] SYM [Al> ... ;Am] (6) 

and is an element of the group algebra s„. 

If we go back to the representation D(S n ), then the action of 5 / fAi,...,A m ] on anv element \v) of V maps this vector 
to an irreducible subrepresentation. A basis of the irreducible representation space can be constructed by applying 
all permutations to 5 / [Ai,...,A m ] \ v ) an d selecting a maximal linearly independent subset. Every standard Young tableau 
gives a different irreducible representation, those corresponding to the same partition [Ai, A m ] but different entries 
are isomorphic. Depending on the structure of D, also the action of the same Young operator on different vectors 
from V can give different irreducible subrepresentations of D. Every irreducible subrepresentation can be constructed 
in the prescribed way. 

Another notion needed in the following is that of the associate representation. For any irreducible representation 
given by a standard Young tableau with partition [Ai, A„] it is given by the transposed standard Young tableau, 
i.e. the rows become the columns and vice versa. The transposed partition is denoted by [Ai, A„] with Ai = m 
and fh = Ai. An example is given in the following figure. 
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FIG. 2. Example for transposing a standard Young tableau 



IV. THE REPLICA SYMMETRIC EIGENSPACES 



We now return to the problem of finding the eigenvalues of the replicated and disorder-averaged transfer matrix 
T n . At the end of the second section we showed that the space Y V® n is invariant with respect to T n for every Y e s„. 



3 



In particular, this is the case for the Young operators which define minimal invariant sets obtainable without further 
knowledge of the exact form of T n . 

In this section we concentrate on a special irreducible subrepresentation described by the standard Young tableau 
with only one row. The Young operator YLi becomes the symmetrizer of the complete symmetric group, its image 
YLiV®" is therefore invariant under permutations. The corresponding irreducible subrepresentations of D are thus 
one-dimensional, all permutations are represented trivially by the identity. Consequently the elements of Y[„] V® n are 
replica symmetric and therefore also the eigenvectors of T n constructed within this space. 

As basis vectors for YLiF®" we introduce 



|Pl, -,pR-l) 



1 



Pi 



! . 



' Pr ! 



E 

=p s Vs=l, 



® \R)® PR 
\s\..s n ) 



(7) 



,R-l} 



where p s , s = 1, R — 1, and pr = n — X^^Li 1 Ps have to be non-negative integers. The replica symmetric submatrix 
of T n can be calculated by 



TtHp 



1, ...,/5r_i |(T1, ...,(Tfl_i) 



(pi, — ,Pfl-i|r n |cri, ...,ctr_i) 



(8) 



(pi, ...,pn_llpl, ...,p R -l) 

The denominator results from the fact that the vectors ([?]) are orthogonal but not normalized. In order to send the 
replica number n to zero we have to introduce generating functions into the eigenvalue equation 



A[„]Z(t7i, (JR-l) = 



E ^(Pl, 

, — ,PR-l} 



by writing 



The eigenvalue equation now reads 



E ''"I 

{pl,...,PJJ-l} 



pi 



.,p/?_i|cri, ...,a R -x)Z{px, 



?-iZ(pi, -,Pr-i)- 



Ar„i$ 



.ri. 



■,XR-1 = 



E 



XcT, 



s-LR,s 



■ $ 



X) s x s rR_i„ 



(9) 



(10) 



(11) 



where we introduced a;^ = 1 for simplicity, for the calculations see appendix A. In this equation a sensible limit n — > 
can be performed. The largest eigenvalue is A = 1 — (3nf + 0{n 2 ) where j3 is the inverse temperature and / the free 
energy. Finally, we change from left to right eigenfunctions, introduce x = (xi, ...,x R -x) G R/^ 1 and 



h r (x) 



ER m 
a=i x S T R,s 

and obtain an equation for a (i? — l)-dimensional invariant density 



$ {0) [x] = J d^y < ^-^(a! - h(y)) > $ (0) [y] 



(12) 



(13) 



where we used the (R — l)-dimensional Dirac distribution <5^ _1 )(-). The density has to be normalized, 
j d(R-i) x $(°)[jc] =10], As in perturbation theory we calculate the 0(n)-corrections of A with the unpertur- 
bated eigenfunction, 



/ = ~ j d^x *<°>[x] « ln(J2x s T R}S ) » 

^ J 8=1 



(14) 



In this paper we do not calculate this free energy for any special distribution. This task itself is very hard and has 
been solved only for a few distributions of quenched disorder, see e.g. and references therein. 

The same equations can be obtained also without using replicas. For the one-dimensional Ising model this was 
established by Derrida and Hilhorst in fllCfl , their method using Riccati variables also generalizes to more complicated 
degrees of freedom than Ising spins. Another result reminiscent of ours was obtained by Lin pjj| , who showed the 
equivalence of an early replica approach by Kac with Dyson's method for the phonon spectrum of a chain of random 
masses and springs. 
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V. TWO-POINT CORRELATIONS 



To be sure that replica symmetry is not violated we have to consider the eigenvalues of the other representations. 
They have a very simple interpretation in terms of connected two-point correlation functions. 
Consider the operator 

X\s) = x s \ S ) (15) 

where x s is any observable assigned to the basis vectors \s), e.g. spin, location, or occupation number. It can be 
simply extended to the replicated vector space V® n by introducing the n operators X^ = l® a_1 gig j®n-a^ 
a = 1, ...,n. They are commutative and measure the value of x at the a-th replica site. In addition we introduce the 
operators 

for any non-negative integer A. For every partition [Ai, A m ] and its transpose [Ai, Am], they can be combined 
to the operator 

X [Al ,...,A m ]=^ (A " l) ®...®^ (X * ) (17) 

acting on V® n . Moreover, it maps any replica symmetric vector to a vector in a representation space belonging to the 
standard Young tableau with m rows of length Ai, A m where we fill one column after the other successively with 
integers 1, n. An example is given by the 2nd tableau in figure 1. A sketch of the proof will be shown in appendix 
B. 

Using this we find that 

tr(T r tX [Al! ... !Am] Tr l X [Al! ... !Am] T^") cx Am] (18) 

for large distances \j — i\. A[ Al A i is the largest eigenvalue if we consider only eigenfunctions in the subspace 
Y[x u ...,x m] V^. """" 

Because of J2 a X a — n we have by definition (|17| ) 

x [Xl ,..., Xm] =x [j:Z2XaM _ Xm] ®i n -T,Z^ . (19) 



Introducing this into (18) we can send n — > and obtain 

^VJ^.a. A..'^--Vv ^ \.. • »>oc l™ A!r^ 2 A.A 2 ,..,A m] > ( 2 °) 

i.e. the two-point correlation function of X^m x x<2 A j decays exponentially with correlation length ^ = 
— 1/lnAj A A2 A ]■ < • > denotes the thermodynamic average in the disordered system with transfer ma- 
trices Ti, i = 1, N. In order to calculate this we still need 2 ^^=2 rea ^ non-interacting replicas of the original 
quenched system. 

Here we concentrate on Young tableaus having only two rows, i.e. to partitions [n — A, A]. There the operator reads 

X [n - X ,x] = (X ® 1 - 1 <8> Xf x <g> 1®(™- 2A ) (21) 

and consequently 

« X [KX] (t) ■ X [KX] (j) » = < (< x Si x Sj > - < x Si >< x Sj >) x >oc A[r A ,' A] (22) 

describes the A-th moment of the connected two-point correlation function with respect to the disorder distribution. 
The correlation length diverges whenever lim„^o Ar n _ A)A i = 1. The criterion for replica symmetry breaking, i.e. the 
degeneracy of the largest replica symmetric eigenvalue with a non-symmetric one, thus coincides with the standard 
criterion for a phase transition. However, for one-dimensional systems with finite R we cannot expect any phase 
transition for non-vanishing temperature. 
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In App. C we will develop an equation for lim„^o Ar„_x,A] ■ The calculations are quite similar to replica symmetric 
one, but due to the more complicated representation structure they are somewhat lengthy. Here we give only the final 
result, an eigenvalue equation for a function 



#[-\,x] . R fi-i ^ R ((*-D A ) (23) 



given by its components ^ 1 X ' X \(x): 



A[-A.Ai*ti A i>(*) - / dR "y £ « s(R - 1] ( x h (y)) II » *iw*(») ( 24 ) 

r 1 ,...,r x = l a=l 

For every eigenfunction of T^ 1-1 ' 1 ' for n — > the function V • is an eigenfunction of the replica 

symmetric transfer matrix given in (|l|) to the same eigenvalue. Only the largest replica symmetric eigenvalue (=1) 
cannot be reached in this way, because the integral of V-$^ _1:1 '(x) over the definition space R fl_1 vanishes. Therefore 
the largest eigenvalue of Tq -1 ' 1 ' equals the second largest of and so on. The first transfer matrix block which 
could produce a diverging correlation length outside the replica symmetric sector is the one corresponding to [—2, 2], 
i.e. to the second moment of the connected two-point correlation function. This is know to be just for spin glass 
transitions where the second moment of the connected two-point function describes the non-linear susceptibility, cf. 

a 

VI. SUMMARY AND OUTLOOK 

In this paper we developed a general replica transfer matrix method capable of handling products of random 
finite-dimensional matrices. We obtained expressions for the free energy (or Lyapunov exponent) from the replica 
symmetric eigenvalues and for the correlation length of the moments of the connected correlation function from non- 
trivial representations of the symmetric group. So we showed that the representation theoretic approach to replica 
symmetry breaking is a general tool for one-dimensional models and probably it can be extended to two-dimensional 
models by considering larger and larger one-dimensional stripes. We have not applied our results to particular 
models, this should done in some future work in order to study the possibility of replica symmetry breaking in one- 
or two-dimensional systems. 

Acknowledgment: I am very grateful to J. Berg and R. Monasson for fruitful discussions and careful reading the 
manuscript. 

APPENDIX A: LAPLACE TRANSFORM OF THE REPLICA SYMMETRIC EIGENVALUE EQUATION 

In this appendix we calculate the Laplace transform of the replica symmetric eigenvalue equation (^|). We start 
with 

A[„]Z(cri, ...,a R -i) = ^2 T^(pi,...,p R - 1 \ai,...,aR-i)Z{pi,...,pR- 1 ) (Al) 

{Pl>— ,PR-l} 

where the replica symmetric transfer matrix is given by 

[Pi, -,PR-1\PU -,PR-l) 

using the replica symmetric vectors 

\pi,...,p R -i)= (A3) 

{s a \J2 l5a«,8=Pa Vs=l,...,iJ-l} 



see Sec. IV. If we introduce the Laplace transformation (|10|) on the left side of (Al) we obtain (introducing xr 
1, pr = n - pi — ... - pr-i) 



G 



A [n] $[x 1 ,...,X R -x]= ^ Z{pi,...,p R -x) ^2 X< 1 ■ -■■ X R R -1 T n l] {Pl,~;PR-lWlT~,VR-l) 
Pl.---.PR-1 gi,---,<tr-i 

= Z(Pi,-,PR-i) ^2^ a5sa - 1 -...-xf^ a - R - 1 (l\^^...^>(R\^T n \s 1 ...s n ) 

pl,...,/P H _l { s °-} 

/ R \ Pl I R \ PR 

E Z (P!> ->PR-l) I ^2 x sTl, s I • - • I ^2x s Tr s J 
px,—,PR-i \s=l / \s=l / 

R \ n R-l /^R rp \P* 

5> r *.. £ % ^ 

s=l / pi,...,p H _i r=l \2^s=l x s±R,sJ 



Pl.---.PR- 

R \ " 



which is equation (Jill). 



> a; s Jfl, s • $ ^ ^ — .•••j-F^ ™ 



APPENDIX B: PROOF TO SEC. V 



(A4) 



In this appendix we show that the operators X[ Xl ,...,x m ] defined in ( p~7| ) map any replica symmetric vector to a 
representation space for an irreducible representation with a Young tableau described by [Ai, A m ]. This can be 
done by proving the equation 

Y[ Xl ,...,x m ] X[x u ...,x m ] F[„] = c [Xl ,...,x m ] X [Xl ,...,\ m ] Y [n] (Bl) 

where C[ Alj Am ] is a real number given by Yj Aj A i = C[ Alj Am ]Y[ Alj Am j. Here we concentrate on the case [n— A, A], 
i.e. to Young tableaus with only two rows. These are the most important cases for our needs, and the proof can be 
generalized directly to more complicated tableaus as well using analogous procedures. 
In the case of two rows we have 

V[n-A,X] = (1 - (1, 2))(1 - (3, 4))...(1 - (2A - 1, 2A)) • SYM [n „ A , A] (B2) 

where (a, b) denotes the transposition permuting a and b, and 

X [n _ KX] = {X <g> 1 - 1 ® X)® x ® i®(»-2A) (B3 ) 

(i) As a first step we note that 

Vtt e S n : tt XW • ... • XW Y [n] = X^, • ... • X^ vr Y [n] 

= 4t)--- X So y H- ( B4 ) 

It follows that SYM[ n _ A A ]A'[ n _ AiA ]Y[ I j] is a sum of certain xl"' • ... • xi"' V'-,, with integer prefactors depending on 
ai < ... < a x . 

(ii) The action of ASYM [n _ A)A] = (1 - (1,2))(1 - (3,4))...(1 - (2A- 1,2A)) on these gives 

ASYM[„_ A)A ] XW ..... X« Y [n] = | y N ^ »p e - 1, M Vp = 1, A (Bg) 

If there were the factors -Y^p-i an d -^p"* f° r an y P — \ the action of (1 — (2p — l,2p)) would annihilate the term. 



The same happens, if there is any p < A for which neither X^lLi nor ^2p^ appear in the product. The sign in (B5) 

(n) (n) 

can be obtained by counting the even indices a p in X„j • ... • JQ A . 

Altogether we find that the action of the Young operator produces only a constant of proportionality, and the proof 
is complete. 
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APPENDIX C: CALCULATION OF NON-SYMMETRIC EIGENVALUE EQUATIONS 



In this appendix we present the calculation of the eigenvalue equations for non-trivial irreducible representations 
at the example of [n— 1,1]. This case is surely the simplest non-trivial one, but the ideas of the calculation are the 
same also for higher representations. 

We consider the standard Young tableau for the partition [n — 1,1] having entries 1, 3, 4, n in the first row and 
2 in the second. The corresponding Young operator 

Y[ n _i il] = (l-(l,2)).SYM(l,3 J 4,...,n) (CI) 

maps an arbitrary basis vector \s ...s n ) up to a normalization constant to 



I* 2 



; o-i, ...,a R -i) := £ (|ss 2 ) - \s 2 s)) ® \ai, <7 mua(S)f , 2) - 1, -, c r maX ( s , s 2) - 1, ctr-i) (C2) 



where the last term in the product is a symmetrized vector in the (n — 2)-fold replicated vector space, cf. equation 
(0), and cr s = X)a=i <V\s- D ue to Y[„_ 11 ]T!, l = Tn^-^i] these vectors form an invariant set with respect to T n . For 
given o"i, o\r_i there are i? — 1 linearly independent vectors of this type, so without loss of generality we can choose 
s 2 = 1,...,R- 1. 

The transfer matrix block to be calculated is 

™[n-i,i]/ i \ (s;cti, ...,cr fl ,_i|T„|r;pi, 

' l (s;f J i: •■•,cr i? _i|r;pi, : 



(s; ai, aR-i\s; at, ctr^i) 
= '5^(TR ! tT 8>r — Tji, r T s ^ t ) 

x ( T »(n-a))[»-2] (£rij ... 5 CTs _ 1; o-^lp!, p r - 1, Pr_x) ■ (03) 

The matrix (t^C™- 2 ))!"- 2 ! j s nothing but the replica symmetric matrix T^"_ 2 2 ' without the average over the quenched 
disorder. For the eigenvalue equation 

A[„_i,i] •C(r;pi,...,p fl _ 1 ) = £ r^ n-1,11 (s;o-i, cr fl ,_i|r; pi, p fl _ 1 )C(s; «7 X) ...,cr fl _i) (C4) 



s;ffi,...,(TH-l 



we introduce again a Laplace transform by 

$ s [x u ...,x R - 1 }= J2 x ° ' ••••••'V • o* i •' • (C5) 

Due to (xr ■= 1) 

£ < : , T^ -1 ' 1 ' (s; (Ti, Cfl-i |r; pi , ar-i) 

Pi,---,Ph-i 

= « £(?Vr s , r - T R , r T s , t ) x f I <' - ■ • - • <~ X ■ - • 4-i (T^"- 2 ))["- 2 l(...|...) J » 

M r \pi,— J 

( R \ n_2 fl-i a; T \ CT9_A 

ty£r \p=l / 9=1 V 2^p=i x p J -R,pJ 

(R \ n R—l 

^x p T R A .^).[]M-r^», (06) 
p=l / r 9=1 

where the 2nd last step is the same as in App. A for the replica symmetric case, and where we use the function h(x) 
defined in (|l2|), the eigenvalue equation becomes 

/ R \ n R-i q, 

&[n-l,l]®r[xi,...,XR-l]= \J2 X P Tr >p) 7T~ $'[ h l(. X )i -> h R -l(x)] . (C7) 

\p=l ) s =l 0xr 
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In the limit n — > this results in equation (24) for A = 1. The calculations for larger A are analogous. 
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